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Abstract 

In this article, the correlated basis-function (CBF) method is applied 
for the first time to the quantum spin-half XY model on the linear chain, 
the square lattice, and the simple cubic lattice. In this treatment of the 
quantum spin-half XY model a Jastrow ansatz is utilised to approximate 
the ground-state wave function. Results for the ground-state energy and 
the sublattice magnetisation are presented, and evidence that the CBF de- 
tects the quantum phase transition point in this model is also presented. 
The CBF results are compared to previous coupled cluster method (CCM) 
results for the spin-half XY model, and the two formalisms are then com- 
pared and contrasted. 

In this paper we consider the T = properties of the quantum spin system 
known as the spin-half XY model, described by the Hamiltonian 
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in the regime < 7 < 1. Note that the index i runs over all N lattice sites 
and that the index j runs over the z nearest-neighbour to i on the linear chain 
(z = 2), the square lattice (z = 4), and the cubic lattice (z = 6). 

In the regime < 7 < 1 the ground state is believed to exhibit Neel ordering 
in the x-direction, and for — 1 < 7 < the ground state is again believed to 
possess Neel ordering in the y-direction. We note that a phase transition point 
occurs for the linear chain model at exactly 7 = 0, and that Neel ordering is 
found to disappear at this point. For spatial dimensionality greater than one, 
the phase transition point of the anisotropic model is also believed to be at (or 
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very near to) 7 = from approximate calculations. The ground state of the 
spin- half XY model on the square and cubic lattices at 7 = is also believed to 
be Neel-ordered in the xy-plane. 

The spin- half XY model was solved exactly by Lieb, Schultz and Mattisjj] 
for the linear chain using the Jordan- Wigner transformation. Since then the 
ground- and excited-state properties have been extensively studied by many 
authors (see, for examples, Refs. [2,3]). However, no exact results exist for 
higher spatial dimensionality, although approximate results such as those from 
spin-wave theory, j§ Monte-Carlo (QMC) methods, §, § series expansions, § 
and the coupled cluster method (CCM) || has proven to be highly successful. 
Extrapolated finite size calculations f| have also been performed for 7 = 0. 

The correlated basis function (CBF) method 10-20 is a widely applied and 
accurate method of modern-day quantum many-body theory. Recently, this 
method has been applied with great success to the Ising model in a transverse 
magnetic field at zero temperature. 21-24 In this article we wish to apply the CBF 
method to the spin-half XY model. We begin this process by firstly performing 
a number of unitary transformations on the local spin axes on two sublattices 
{A, B} in order to simplify the problem. The first such transformation on the 
A-sublattice is given by, 

a x -> a z ■ a y -> a x ; a z -f a y , (2) 

and the second transformation on the _B-sublattice is given by, 

a — ► — it ; <j a — > —a ; a — > a" . yd) 

Note that both of these transformations are simply rotations of the local spin- 
axes of the spins, and that the eigenvalue spectrum of the problem is left un- 
changed because these transformations are unitary. The Hamiltonian may now 
be rewritten in terms of these new spin-axes as 



H 

8 



(i-7KX + (i + 7 K^ 



(4) 



We may now define a ground-state trial wave function, given by 

1 N 

\i/>) = cxp{[/} |0> ; tf=25>(r«)<7fo? , (5) 

i<j 

where u(iY, ) is the pseudopotential. The reference state |0) is given by a tensor 
product of spin states which have eigenvalues of +1 with respect to tr z , and this 
state is an exact ground eigenstate of the Hamiltonian Eq. (|j) when 7 = 1. 
Translational invariance also implies that the pseudopotential, u(ry), depends 
only on the relative distance, n = — r, = r^. 
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The treatment of the spin-half XY model by the CBF method is continued 
by defining the lattice magnetisation (i.e., again the magnetisation in the z- 
direction in terms of the rotated local spin-axes), given by 

(ip I of I ib) 

for a ground-state trial wave function, \tp). Furthermore, the 'transverse' mag- 
netisation (in terms of the rotated local spin-axes) is given by, 

_ {j, | of | jj) 

We may now define a spatial distribution function (which plays a crucial part 
in any CBF calculation) in the following manner, 

(ib I erf a 3 ; I ib) 

Furthermore, we may also determine an expression for the expectation value of 
the ground-state energy of the spin-half XY Hamiltonian of Eq. (Q), where 

e = (mjj) 

N N{ib~\ib~) ' [ ' 

In the region 7 > we now make the explicit assumption that A — 0, which 
is in agreement with our Ansatz for the trial wave function of Eq. (|J). This 
assumption furthermore implies that A/N = (P - *) — (P*~) = — = 0, 
where P^ and P* - are spin projection operators in the positive and negative 
re-directions respectively. However, we note that + = 1 must also be 
correct, which therefore implies that p^ = p^~ = 1/2. Hence, we may treat this 
problem completely analogously to a binary-mixture of two types of bosons [p5[ 
each with a density equal to one-half. 

The expression in Eq. (^) may be determined via a hyper-netted chain 
(HNC) cluster expansion, and it is readily found using this procedure that the 
ground-state energy is given in terms of a functional with respect to the pseu- 
dopotential, it(n), where 
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Note that A(n) is unity if n is a nearest-neighbour vector and is zero otherwise. 

Self-consistent HNC equations may also be determined. These equations 
may be then iteratively solved, and thus G(n) (and so the ground-state energy) 
may be also obtained. The first method of finding the pseudopotential has a 
"variational" flavour, and we parametrise it(n) in the following way, 

u(n) = a A(n) . (11) 
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Figure 1: Results for the ground-state energy of the spin-half XY model for 
the square lattice plotted as a function of the strength of the nearest-neighbour 
pseudopotential, a, for varying 7. At 7 = —0.36 we see that the minimal 
solution that we have tracked from 7 = is lost. 



A(n) is unity if n is a nearest-neighbour vector and is zero otherwise. We 
now minimise the ground-state energy with respect to a at a given value of 7. 
Indeed, at 7 = 1 we already know that all correlations have zero strength as 
our reference state |0) is an exact ground eigenstate of Eq. (||), and this implies 
that a = 0. We thus track this solution at 7 = 1 in the regime 7 < 1, and the 
ground-state energy as a function of a for various values of 7 is plotted in Fig. 
for the square lattice. We may see that at 7 = —0.36 the minima that we 
have tracked from 7=1 become a point of inflection. 

The second such method of determining the pseudopotential is to determine 
the optimal value for the function u(n) with respect to the ground-state energy, 
E/N . This is stated as, 

ou(n) 

which may be determined analytically from Eq. ([l0|). In the context of this 
article, this approach shall be referred to as the P aired- Phonon Approximation 
(PPA), in analogy with a binary mixture of two types of bosons, for example. 
Note that we do not explicitly state here the resulting PPA equations for this 
model, although the treatment is fully analogous to that performed for the 
transverse Ising model and the interested reader is referred to Ref. [22] for a 
full account of this calculation. 

For details of the specific application of the CCM to the spin-half XY model 
the interested reader is referred to Ref. [8]. We note however that two types 
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Figure 2: CBF results for the ground-state energy of the spin-half XY model 
on the square lattice compared to results of high-order CCM results of Rcf . [8] . 

of approximations are made, namely, the SUB2 approximation which retains 
all two-body correlations in the approximate CCM ground-state wave function, 
and the LSUBm which retains all correlations in a locale defined by m. 

Results for the CBF ground-state energy of the spin-half XY model on the 
square lattice compared to results of high-order CCM results are given in Fig. 
H and, for the isotropic point (7 = 0) only, in Table ||. We may see from Fig. || 
that both sets of results are in excellent qualitative agreement over a wide range 
of 7. It is furthermore seen from Tables that CBF results are in excellent 
quantitative agreement with LSUB2 CCM results at 7 = 0. This is a perfectly 
reasonable result because both the CBF and CCM LSUB2 results only utilise 
two-body correlations. It is, however, expected that the inclusion of higher- 
order correlations in the CBF trial wave function would produce more accurate 
results for the energy, as is seen for the CCM. Thus, from Tables |-|, we see 
that the CBF results at 7 = capture about 59% of the correlation energy 
for the linear chain, 76% of the correlation energy for the square lattice and 
85% of the correlation energy for the cubic lattice (in comparison with exact 
and extrapolated CCM results). Indeed, the extrapolated CCM results present 
some of the most accurate results yet seen for the isotropic XY model on the 
square and cubic lattices. (Results for the linear chain and cubic lattice are 
qualitatively similar to the results presented for the square lattice in Fig. || and 
so are not plotted here.) 

Results for the sublattice magnetisation of the spin-half XY model on the 
square lattice are presented in Fig. || and in Table || for the isotropic point, 
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Figure 3: CBF results for the sublattice magnetisation of the spin-half XFmodel 
on the square lattice compared to results of high-order CCM results of Ref. [8] . 



7 = 0. Again, it is seen from Fig. g that the CBF results are in good qualitative 
agreement with the known results of this model. However, the CBF result for 
the sublattice magnetisation at the isotropic point (7 = 0) is slightly too high, 
although it is again expected that higher accuracy would be achieved with the 
inclusion of higher-order correlations in the approximate CBF ground-state wave 
function. Again, results for the linear chain and cubic lattice are fully analogous 
to the square lattice case and so are presented only for the isotropic model in 
Tables | and |. 

Results for the phase transitions points predicted by the CCM method are 
also given in Tables although no such results are explicitly given for the 
CBF method in these tables. It is however noted here that the loss of "minima" 
within the parametrized HNC CBF approach (at 7 = —0.36 for the square 
lattice) may be associated with a phase transition within this system. This 
constitutes a powerful result for such a simple variational-style calculation. Note 
that similar behaviour is also seen for both the linear chain and cubic lattices. 
An analogous change in the energy surface with respect to u(n) for the CBF 
PPA approach seems to occur for varying values of 7. However, in this case, 
the situation is much less clear-cut because, near to this point, convergence of 
the PPA equations becomes very difficult. 

In this article, the CBF method has been applied with much success to the 
quantum spin-half XY model on the linear chain, the square lattice, and the cu- 
bic lattice in order to obtain accurate results for the ground-state energy and the 
sublattice magnetisation. These results were found to be in excellent qualita- 
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Table 1: Ground-state energy and sublattice magnetisation for the one- 
dimensional XY model at 7 = compared to exact results of Ref. [3] and 
CCM results of Ref. [8] . The critical values of 7 for the anisotropic model are 
also given. 
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tive agreement with previous CCM calculations, |8J although more quantitatively 
accurate CBF results would be possible with the inclusion of higher-order cor- 
relations (than Jastrow correlations) in the trial ground-state wave function. A 
strength of the CBF method is that it is not limited by the presence of frustra- 
tion, in contrast with QMC methods for example, and a direct extension of this 
work would be to include next-neighbour-neighbour interactions in our model. 
Indeed, this presents the possibility that this calculation might be utilised to 
provide a trial or guiding wave function for these QMC techniques in the pres- 
ence of such frustrating next-nearest-neighbour bonds. Also, it is possible to see 
that an extension of this work to Heisenberg antiferromagnetic (HAF) models 
could follow a similar path to that outlined in this article. One would per- 
form a similar set of rotations of the local spin-axes and then perform HNC 
re-summations of the relevant quantities that one is interested in. For exam- 
ple, one might consider the HAF on the triangular lattice, and in this case one 
would perform a rotation on three sublattices - such as that utilised by Singh 
and Huse for this model. |27j] Previous CCM results for the spin- half XY model 
quoted in this article were also seen to provide excellent results for this model, 
and they are furthermore a valuable yardstick with which to compare our new 
CBF results with. 

We note that the CBF approach utilises a Jastrow wave function and its 
bra states are always the explicit Hermitian adjoint of the corresponding ket 
state. Hence, for the CBF approach, an upper bound to the true ground-state 
energy is, in principle, obtainable, although the approximations made in calcu- 
lating the energy may destroy it. By contrast, the CCM uses a bi- variational 
approach in which the bra and ket states are not manifestly constrained to be 
Hermitian adjoints and hence an upper bound to the true ground-state energy is 
not necessarily obtained. Also, the CCM uses creation operators with respect to 
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Table 2: Ground-state energy and sublattice magnetisation for the square lattice 
XY model at 7 = compared to CCM calculations of Rcf. [8] and series 
expansion calculations of Ref. [7]. The critical values of 7 for the anisotropic 
model are also given, where the value in parentheses is the estimated error in 
the final decimal place shown. 
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some suitably normalised model state in order to span the complete set of (here) 
Ising states. The CBF method, in essence, uses projection operators to form 
the Jastrow correlations with respect to a reference state, |0). In some sense, 
the CCM is found to contain less correlations than the others at 'equivalent' 
levels of approximation (e.g., the CCM LSUB2 approximation versus Hartree 
and nearest-neighbour Jastrow correlations). A strength of the CCM is that 
it is well-suited to the inclusion of high-order correlations in the approximate 
ground-state wavefunction (for example, via computational techniques). Fur- 
thermore, the CCM requires no information other than the approximation in 
S and S in order to determine an approximate ground state of a given system. 
The CBF method, however, may require that only a certain subset of all possible 
diagrams are summed over (e.g., the HNC/0 approximation). 

We finally note that the results of one method reinforce and sometimes 
elucidate the results of the other, and the application, in parallel, of two such 
methods to the same model can lead to a deeper understanding of the behaviour 
of it. 
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